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We discuss the possible existence a spin singlet excitation with charge ±2 (p-mode) originating the Aig 
Raman resonance in cuprate superconductors. This Tj-mode relates the c/-wave superconducting singlet pairing 
channel to a d-wave charge channel. We show that the rj boson forms a particle-particle bound state below the 
2A threshold of the particle-hole continuum where A is the maximum d-wave gap. Within a generalized random 
phase approximation and Bethe-Salpether approximation study, we find that this mode has energies similar to 
the resonance observed by Inelastic Neutron Scattering (INS) below the superconducting (SC) coherent peak at 
2A in various SC cuprates compounds. We show that it is a very good candidate for the resonance observed in 
Raman scattering below the 2A peak in the Aig symmetry. Since the T]-mode sits in the 5 = 0 channel, it may 
be observable via Raman, X -ray or Electron Energy Loss Spectroscopy probes. 


I. INTRODUCTION 

The study of collective excitations of a condensed matter 
system can provide informations about the symmetries of an 
order parameter or the nature of the interactions between par¬ 
ticles. The analysis of collective modes can also help us to 
understand the origin of high critical temperature supercon¬ 
ductivity (SC) in cuprate compounds (see schematic cuprate 
compounds phase diagram on Fig. 1). 

The existence of a collective spin excitation in the SC state 
of cuprates has been highlighted by Inelastic Neutron scat¬ 
tering (INS) experiments with the observation of a resonance 
at 41meV in YBCO compounds and at similar energies in 
other compounds*"^ around the antiferTomagnetic (AF) order¬ 
ing vector Q= {7t,n). This resonance, which stands below 
the threshold of the particle-hole continuum (A^^^ is the 
maximum of the ^f-wave SC gap), seems to scale with the su¬ 
perconducting gap energy above optimal doping (at least until 
p = 0.19 hole doping in Bi2212)^. 

This INS resonance has initially been explained in the 
framework of the SO(5) emergent symmetry model for 
cuprates as a ;r-mode^’®. However this ;r-mode lies at higher 
energies than experimentally observed, and the neutron reso¬ 
nance is now explained as a spin triplet exciton or resonance 
which emerges in the SC state because of a residual spin-spin 
interaction in the system^ (see also Ref.^). 

A resonance very similar to the Neutron resonance has also 
been observed in Raman scattering experiments in the Aig 
symmetry channel"^"*^. In YBCO the Ajg Raman resonance 
is located at 41 meV at optimal doping, and follows the neu¬ 
tron resonance energy with nickel and zinc substitutions*^’*^. 
This resonance is not seen in the Big symmetry, which is scan¬ 
ning the anti-nodal region of the Fermi surface, nor in the B 2 g 
channel, which is scanning the nodal region. In the Big chan¬ 
nel, a SC coherence peak is observed at higher energy than the 
Aig resonance. Its energy matches well twice the maximum 
of the d-wave SC gap 2Aj(^ observed in other spectroscopies. 

The Aig Raman resonance has initially been attributed to 


the 2A^p SC pairing gap. However, later considerations 
showed that long range Coulomb screening washes out the 
single particle pair-breaking contribution to the Raman re¬ 
sponse in the Aig channel leaving the position and the inten¬ 
sity of the Aig Raman resonance essentially unexplained*®"*®. 
Note that the two-magnons process proposed in Ref.*^ pro¬ 
duces a resonance at twice the energy of the neutron resonance 
and thus cannot explain the near perfect energy matching be¬ 
tween the Aig Raman resonance and the Neutron resonance 
around the optimal doping. 

From a theoretical point of view there are several open 
questions regarding the Raman resonance in the Aig channel: 
Since the Aig symmetry scans the whole Brillouin Zone (BZ), 
how can we explain the absence of a superconducting coher¬ 
ent peak in the Aig channel whereas it is observed in the Big 
channel? How can we explain that the Aig resonance emerges 
at lower energy than the SC coherent peak one? How can we 
explain the perfect match between the energy of the Raman 
resonance in the Aig channel and the INS frequency resonance 
at optimal doping? What are the informations that this Aig Ra¬ 
man resonance could give us about the physics of Cuprates ? 

In this paper, we present a coherent scenario that provides 
an answer to the latter questions. Our scenario builds on re¬ 
cent experimental developments^**"^** that pinpointed the ex¬ 
istence of a Charge Density Wave (CDW) state, with order¬ 
ing vector qc= {qx,0) and qc= {0,qy) with qx = qy « 0.3;r 
which exists between hole doping p « 0.1 until p « 0.19^* 
(see Fig. 1, note that the endpoint of the CDW order is still un¬ 
der debate). We explain the Raman resonance in the Aig sym¬ 
metry as a collective mode that allows excitations between 
the £/-wave SC pairing sector and the observed tZ-wave CDW 
charge sector. The proposed scenario is supported by the fact 
that the Aig Raman resonance exists between p — 0.12 and 
p — 0.22^^ which overlaps with the hole doping range where 
a coexistence between a rZ-wave SC and a tZ-wave charge den¬ 
sity wave (CDW) ordered states have been observed^* (see 
figure 1). This collective mode is a spin singlet {S = 0) exci¬ 
tation with a charge ±2 and we call it a Tj-mode of the system. 
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FIG. 1. (Color online) Schematic Temperature-hole doping (T,p) 
phase diagram of superconducting cuprates. Close to zero hole dop¬ 
ing, the system exhibits a Mott insulator antiferromagnetic phase 
(AF). At higher hole doping, the Pseudo-Gap (PG) phase appears 
below the critical temperature T* (dashed line). A dome of a d- 
wave superconducting (SC) phase appears at average doping below 
the critical temperature Tc (dotted line). The maximum of critical 
temperature Tc is obtained at the optimal doping pg. Above T* at av¬ 
erage doping, the cuprate compounds exhibits a metallic phase with a 
linear temperature dependence of the resistivity called Strange metal 
(SM). The conventional Fermi liquid (FL) state appears at high dop¬ 
ing. A charge density wave (CDW) order appears below T* and 
above Tg (blue area). The hole doping where the CDW state disap¬ 
pear is very discussed then we do not finish the transition line. Note 
that T* meet Tg around p = 0.22 in Bi2212 compound only. In other 
compounds, the behaviour of the T* line can be different. 


In order to explain the Raman resonance, we assume the 
coexistence between a tZ-wave SC phase and a tZ-wave CDW 
ordered phase. The CDW phase has the same ordering vector 
as the observed one qc= {qx,0) and qc= {0,qy) with qx = 
qy « Q3n. This coexisting phase has been observed by STM 
in Bi2212 compound^*. 

We note that in the references^^’^"*, the authors develops an 
effective t — J model with coexistent SC and CDW order that 
reproduce qualitatively the dependence of the frequency res¬ 
onance in the underdoped phase. However, they considered 
a charge ordering wave-vector Q = {7t,7z), which is not ob¬ 
served experimentally and did not reproduce the peak in the 
Big symmetry. 

Here, we argue that the Ai^, Raman resonance is the collec¬ 
tive mode describing the proximity of the charge and pairing 
channels. Fluctuations between those two sectors are typically 
described by non-linear a- models as in Refs^^’^®. Within a 
model of itinerant electrons interacting through an effective 
AF spin-spin coupling, we find that the Tj-resonance forms a 
particle-particle bound state situated below the 2Ao-threshold, 
in a very similar way to the triplet spin exciton revealed by 
neutron scattering. 

The paper is organized as follows: in the section II, we 
present the theoretical model we use to reproduce the experi¬ 
mental Raman data. We explain how we describe the different 
phases of the Cuprates (section II A). We present the differ¬ 
ent susceptibilities we consider (section IIB) and give details 


about the evaluation of the different Feynman diagrams (sec¬ 
tion IIC). We also explain the calculations of the full Raman 
response (section IID). In the section III, we present the re¬ 
sults of our theory and we compare them to the experimental 
data before to conclude (section IV). Explanations and infor¬ 
mation about the spin fermion model and the experimental 
setup are given in details in the appendix A and B respectively. 


II. THEORETICAL MODEL 

A global explanation of the phase diagram of cuprates com¬ 
pounds is still lacking. For example, we still do not know 
whether the cuprates are fundamentally doped Mott insulators 
admitting a Coulomb energy repulsion of leV which is crucial 
to explain the emergence of superconductivity (see e.g'^^“^®), 
or whether an itinerant electron picture with strong Anti- 
Ferromagnetic (AF) fluctuations is a good approximation to 
explain the main features of the phase diagram (e.g.^*^^). A 
consensus exists in the recognition that three main players are 
present in the phase diagram: tZ-wave superconductivity, AF 
order and fluctuations, and Mott insulating phase. Note that 
the charge order could be a fourth key player of the cuprate 
compounds physics. 

Below T*, the PG phase critical temperature'*^"^^ (see Fig. 
1), the influence of charge orderings seems to be stronger 
than previously thought. The proximity between pairing and 
charge channel is well described by spin fermion approach. 
Such approach have been used through simplified Eight Hot 
Spots (EHS) models in order to explain the PG phase^®’^®. In 
Ref a coexisting d-wave bond order and li-wave SC state 
are related by SU(2) symmetry. Note that the SU(2) symmetry 
relating the charge and spin channels is only realized in a strict 
sense in the eight hot spots model where the Eermi dispersion 
has been linearized around each hot spot. 

However, the CDW and PG phases are different phases as 
the CDW phase appears on the tips of the Eermi arcs after the 
antinodal zone has been gaped out (see Eig. 1). The study 
of SC and CDW competing orders related by angular fluctua¬ 
tions was performed in a model based on Landau theory^^. In 
this model, the SU(2) symmetry can be recovered by ignoring 
the quartic and anisotropic couplings. Note that in Ref.^^, a 
spin fermion approach was developped to explain the stabil¬ 
ity of CDW phase related to modulated SC order called pair 
density wave (PDW) which has not been observed experimen¬ 
tally. In the following, we leave aside the complexity of the 
spin fermion approach and we describe the SC/CDW coexist¬ 
ing phases with another model. 


A. Description of the system 

1. The two body Hamiltonian 

We consider a system of itinerant fermions interacting 
through an effective AE spin-spin coupling close to vector 
Q derived from spin-fermion approach^® (detailed in Ap¬ 
pendix A). We then describe a system of fermions with a ki- 






3 


netic energy and nearest-neighbor (n.n.) interaction Hamil¬ 
tonian involving both spin and charge in the weak cou¬ 
pling regime. The Hamiltonian writes H = Y.i.j.atijCjfjCja + 
Y.(^ij){JijSi-Sj + Vijninj) where n; = and S,- = 

are the density and spin operators respectively 
with Gap the Pauli matrix vector, {i j) denotes summation over 
nearest neighbors and f,y is the hopping parameter. Jij is the 
n. n. super exchange coupling and Vjj the n.n. Coulomb term. 
Long-range Coulomb effects will be considered later while 
discussing the Raman response. We will neglect them for the 
study of the collective mode. Involving the Fourier transform 
Ci.a = where N is the total number of lattice 

sites, the Hamiltonian writes; 


k,(7 

+ Yi (■^q‘'k.a‘^a/3‘'k+q,/tCk'+q,y‘^}'5‘^k',5) 
k.k'.q 

+ Yj ^q‘'k.<T‘'k+q.<rCk'.(T'‘'k'-q.o''; 

k.k'.a.fj' 


where Cp^a is the annihilation (creation) operator of an elec¬ 
tron with spin a and impulsion p. is the one particle en¬ 
ergy with a tight-binding form obtained to fit ARPES data'*^ 
: i^k = —2t{cos{kxa) -t-cos(kya)) + 4t'cos{kxa)cos{kya) + 
to{cos{kxa) — cos{kya))^ — fX where t,t' are respectively the 
first and second neighbor hopping terms with t' = —0.3f 
and to — 0.084f which is third and fourth neighbor hopping 
term, a is the cell parameter set to unity and ji is the chem¬ 
ical potential determined to adjust the hole doping. Jq — 
— 5 {cos{qxa) + cos{qya)) is the Fourier transform of Jij de¬ 
veloped around Q, while Vq = j {cos{qxa) + cos{qya)) de¬ 
notes the amplitude of the nearest-neighbour Coulomb inter¬ 
action for small values of q. Vq accounts for repulsion be¬ 
tween charge and vanishes in conventional t — J model. In the 
following, we treat Vq and Jq as independent parameters'’'^®. 


2. Order parameters and ejfective Hamiltonian 


In the following, we focus on the superconducting phase 
where the Raman Aresonance peak has been observed. We 
consider that the superconducting phase around optimal dop¬ 
ing is a coexistence between a rZ-wave SC order and the d- 
wave CDW order. Thus, we decouple the Hamiltonian (1) in 
the SC pairing channel and in the charge channel by intro¬ 
ducing the two order parameters Ajc.k and Acdwm- ^ 5 C,k is 
the superconducting order parameter which describes the d- 
wave SC order. Acow.k is the CDW order parameter which 
describes the tZ-wave charge order with ordering vector Qp. 
Applying the Hubbard-Stratanovitch transformation to the 
Hamiltonian (1), the effective Hamiltonian of the system can 

be written in the basis = ( 4 ,c,,c_k-qc.cf, 4 +qc.cy>‘^-k.cf) 


as 


(k,®) 


/ Zcu-^k 
0 
At 

“CW,k 
\ ^5C.k 


0 

T “s-k-qc 
t^5C,k+qc 

_ At 

^CDW.k 


^CDW,k 

At 

^5C.k+qc 
~ ^k+qc 
0 


^SC.k \ 
— ^CDW.k 
0 

ico + ^-k ) 
( 2 ) 


with ^_k = ^k- Note that the Green function inverse matrix 
is very similar to the one calculated in^^’^'^. From invert¬ 
ing the matrix (2) , one can find the Green functions of the 
system. Particularly, the CDW quasiparticle Green function 
(jAcw (jj and the anomal SC Green function (jj) 

writes : 


G^®u/(k^ej)=Acw,kX 

(^L,k+qc^SC,k + |AcW,kP - ^k+qc^t: + 0)^ - Za)(i*k+qe + “^k)) 
det (G(k,a))) ” 

(3) 




k+qc 


|Asc.k+qc 1 ^] Asc,k + Ajc.k+qc |Ac£IW,k 


det(G(k,a))) 


(4) 


with 

det(G^^ (k, ®)) = |AcDW,k|^ 

+ 2 |AczJW.kP (A^c.kAsc.k+qc + — ^k+q^^k) 

+ + ^k + |Asc,kP) + “^k+qc + I Asc,k+qc | ) • (5) 

The Green functions of the system will be crucial to estimate 
the different susceptibilities of the system (see section II B). 


3. Phenomenology and symmetries of the SC and CDW orders 

In principle, we have to solve the self-consistent equations 
to determine Asc{k,Co) and AcDw{k,Co)- However, the goal 
of this paper is not to explain the (T,p) phase diagram. Based 
on experimental evidence that the charge is modulated inside 
the superconducting phase of cuprates^*,we assume that the 
SC phase is composed by a competition between a cZ-wave SC 
phase and a cZ-wave CDW order. The amplitude of the order 
parameters A^c.k and AcDW,k is expected to vary with the hole 
doping. The CDW order has the same ordering wave-vector 
Qc as the one observed experimentally. In this paper, we do 
not address the problem of the PG state. We then neglect the 
PG energy scale and we assume that the PG state does not 
infer with the Raman response in the SC state. This assump¬ 
tion is expected to be valid around optimal doping and in the 
overdoped regime. 

In Bi2212, the SC order parameter decreases above optimal 
doping p = 0.16 and vanishes around doping p = 0.25 as men¬ 
tioned in^^ The CDW order is maximal at doping p = 0.12 
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and vanishes around p = 0.19 at zero temperature^^ while the 
PG transition line meet the superconducting transition line at 
a slightly higher doping, around p = 0 . 22 ^''’^^ 

We assume a simple momentum and fre¬ 
quency dependence of the SC pairing Asc{k,(i>) = 

{cos{kxa) — cos{kya)) f{co,A^Q) and CDW order 
AcDw{k,co) = ^^{cos{kxa)-cos{kya))f{co,A‘^iy^) 



_ 

order and f[( 0 ,X) = e li? with a variance CJ = 2f/1.177 
ensuring a half width at half maximum equals to X. Recent 
STM experiments demonstrate the cZ-wave character of CDW 
order^^. 

Above optimal doping, the phase diagram (fig 1) shows that 
the superconducting state disappears at an higher temperature 
than the CDW state and we can consider that the CDW gap 
amplitude is lower than the SC one : In the 

underdoped regime, the decreasing of Tc which goes below 
the CDW critical temperature (see fig 1) can be described by 
the opposite regime A^^^ < A®^^. 

The frequency dependence of the order parameters is a con¬ 
sequence of the spin-fermion model as studied in ref.^®’^^ and 
in Appendix A. We choose a frequency dependence of the or¬ 
der parameters that is a good approximation of the one coming 
from the spin fermion model (see Appendix A). 


4. SU(2) symmetry between SC and CDW order 

The concept of emerging symmetries has first been intro¬ 
duced to explain PG phase of cuprates compounds. Particu¬ 
larly, in Ref.^^’^^, a representation with a pseudo-spin opera¬ 
tors was introduced which rotates the cZ-wave SC state onto 
a charge order with ordering vector qc = In systems 

with coexisting ^Z-wave SC and CDW states, one can define 
a pseudo-spin respecting the SU(2) symmetry which rotates 
from SC state onto a CDW order^®. Indeed, in d-vtave SC and 
CDW coexisting orders, we can define a pseudo-spin where 
each components is a SC or CDW operator with Aq = Acdw, 
Ai = A^^and A_i = — A|, 


(6a) 

■^0 = X 52 ^kCkfj'-k+qcO') (6b) 

^ k.fj 

A_i = ^^4Ck;C_k|. (6c) 

v2 k 

The lowering and raising pseudo-spin operators 77 +, t]^ = 
(77 + )^, and 77 ^, follow the definition ; 

^ “ 524t‘'-k+qci 

k 

= 52 (“^kt^kt + ‘^k+qci'^k+qci “O’ C^^) 

k ^ ^ 


The operators (7 form an SU(2) Lie Algebra. The SU(2) spin 
operators (7) rotate each component of the pseudo-spin ( 6 ) 
into another^®. 

Typical approach to study this underlying SU(2) sym¬ 
metry are non-linear a models^^’^®. Such models 
have been proposed to explain the physics of cuprate 
superconductors^^^’^®’^®’^^. In the following, we calculate the 
response of the operators (7a) that connect SC and CDW or¬ 
ders. 


B. Susceptibilities of the system 

J. The collective modes 

In order to study the collective mode physics, we analyze 
the linear response for the two following operators. First, the 
INS resonance originates the excitation from the SC singlet 
state to the spin triplet at wave vector Q associated with the 
AF order^. It is described by the operator 

5+=A-'/2^4,^Ck+Q;. (8) 

k 

which destroy a bosonic excitation at momentum Q with a 0 
charge and spin 1 as presented in Fig. 2 d). The 77 mode 
excites from the SC order parameter (Ckt‘'2k|) '■ 6 ® CDW 

order (Ck^Ck+qci). It is described by the operator: 

77 =7V-'/2£ck+q,tC-ki, (9) 

k 

where qc is the ordering vector of the charge density wave 
(CDW) order. The 77 operator destroy a bosonic excitation at 
momentum Qc with a 0 spin, a —2 charge and is fully symmet¬ 
ric as presented in Fig. 2 e). Upon action of it, the SC state 
transforms into the CDW state in the same way than the op¬ 
erators (7a). The 77 -mode is obtained at wave vector Qc with 
Qc = (0, qy) = {qx, 0) and qx = qy = 03%, contrasting with the 
typical Q = {%, %) location of the spin triplet exciton. Within 
the framework of the eight hot spots spin fermion model, 
where the electronic dispersion has been linearized around the 
Fermi points, the operator described in Eqn.(9) represents an 
SU(2) rotation between the CDW and SC states. 


2. Susceptibilities in the Random Phase Approximation ( RPA) and 
Bethe-Salpether Approximation 

We rewrite the Hamiltonian (1) with the spin operators (S± 
of equation (8)) and the 77 -mode operator (see equation (9)) 
: H = Lk.CT^kCkcrCka +-/55+5_ +Jrirj^rj. The magnitude of 
the interactions are 7^ —3J — V in the 77 -mode channel and 
Js = 2J in the spin channel. We compare the spin suscepti¬ 
bility derived from the Random Phase Approximation (RPA) 
with Xs = —id (0 ([‘5+ (0 (0)]) depicted in Fig 3 and 

the 77 -mode susceptibility derived from Bethe-Salpeter ladder 
approximation with x-q = —id (f) ([t? (0 i (0)]) shown in 
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k . kk + Oc . k—k 

I J (Jfnorm — _^_ ^^CDW = (^^sc — ^ O ^ 

FIG. 2. (Color online) Graphical representation for a) the photon- 
electron vertex at q and for the symmetry A and b) the neutron-spin 
vertex at Q = (k, k). Graphical representation for a) the Coulombian 
interaction vertex (see ref.'^ for details), b) the spin triplet exciton 
vertex and c) the T]-mode vertex. In f) are presented the diagram¬ 
matic representation of the normal Green functions, the CDW Green 
function and the SC anomal Green function. 


a) ^ 

k 



k + q + qc k + q + q„ + qi k + q + qe + qi+qz 

- ■ - > -[□->— 

+ "F... 

— ■ - > -□ —>— 

— k -k - qi -k - qi - q2 

FIG. 4. (color online) a) The Feynman diagram of the r\ collective 
mode. The red circle represents the CDW quasiparticle Green func¬ 
tions of equation (3) . The blue circle represents the anomal 

SC function of equation (4). is the Raman vertex in the A 
symmetry, b) T is the diagram resulting from the Bethe-Salpether 
series of diagram. £„qn is the summation over the internal momenta 
qn. qc is the charge density ware ordering vector. The value of the 
Tj-collective mode susceptibility is evaluated by the formula (10) 
in the approximation of a frequency and momentum independent . 



FIG. 3. (color online) a) The bare polarization bubble we use for 
the spin mode. For the spin mode, the external momentum equals the 
AF ordering wave vector Q = (tt, k). b) Series of Feynman diagram 
of the spin mode resulting from the RPA approximation. The spin 
susceptibility Xs is the result of this series of diagram which can be 
evaluated by the formula 10. 


Fig 4. Within the RPA and Bethe Salpether approximation, 
the fall susceptibilities are given by; 




Xn\co) 


Xs{(o) 


l-F75A5(cd)’ 


( 10 ) 


where Xg is the bare polarization bubble for the spin mode (see 
Fig. 3 (a)), Xri is the single segment of the ladder diagrams 

describing the Tj-mode (see Fig. 5 (a)) and Xri is the hrst 
term of the diagram series for Tj-mode (see Fig. 5 (b)) in the 
A symmetry. Note that the detailed calculation of the bare 
polarization bubble and diagrams is done in the section IIC. 


C. Evaluation of the diagrams 

In this section, we focus on the diagram calculation pre¬ 
sented in Fig. 3 (a), for the bare polarization of the spin mode 


k + q + qc k + q + qc+ qi 

—9— 



FIG. 5. (color online) a) Xrj is the single segment of the ladder dia¬ 
grams describing the 77 -mode. This term is evaluated in the formula 
(12) in the section IIC 1. b) Xq^’^ is the first term of the serie of 
ladder diagram, is the Raman vertex in the A symmetry. This 
term is evaluated in the formula (14) in the section IIC 2. Both terms 
are evaluated in the approximation of a momentum and frequency 
independent Jrf. 


a) X° = 


Xs, in Fig. 5 (a), for the single segment of the ladder diagrams 
describing theti-mode x^ and in Fig. 5 (b) for the first term of 
the T 7 -mode diagram series. The explicit internal and external 
momenta dependence of each diagram are explicitly presented 
in Fig. 3 (a) and 5. 

In order to include the frequency dependence of the or¬ 
der parameters in the diagram calculation, we develop the or¬ 
der parameters depending on the internal frequency standing 
inside the summation terms around the external frequency, 
A(e,„,) « A(e„f) + (e,>„ - and we keep 

only the zeroth order term (with eint(ext) the internal (exter¬ 
nal) Matsubara frequency). This implies that our susceptibil¬ 
ity will be under-estimated at high-frequency. 
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SC/CDW 



- 0.2 - 0.1 0.0 0.1 0.2 

co(eV) 


FIG. 6 . (Color online) Imaginary part (dashed line) and real part 
(solid line) of (a) the spin bare function at q = Q = (k, k) and 
(b) the 77 -mode bare function Xri ^t q =(0,0) atT = OK in a d-wave 
superconductor with = 60meV mixed to a d-wave charge order 
with = 35meV . A broadening 5 of 5meV was employed. 


We also assume that the T 7 -mode interaction Jri and the spin 
interaction Js do not depend on the frequency and the momen¬ 
tum. This approximation allow us to decompose the ladder 
diagram of Fig. 4 in the diagram of Fig. 5. The momen¬ 
tum sum is done over the first Brillouin zone with meshes of 
400 X 400 after doing the summation over internal Matsubara 
frequencies at zero temperature T = 0 . We do the analyti¬ 
cal continuation on the external Matsubara frequency writing 
ico = (0 + i5 in the energy denominators. A small broadening 
is introduced by the parameter 5 and can be understood as a 
residual scattering. 


1. Calculation of the bare polarization bubbles. 

The bare polarization is the response function of non lo¬ 
cal density operator : = ~'®(F)(Pq(F)pq( 0 )) with Pq = 

l/^LkCfcCk+q. From the figure 3 (a) and following the aca¬ 
demical way to evaluate these diagrams (see for example the 


Fermi Liquid 



FIG. 7. (Color online) Imaginary part (dashed line) and real part 
(solid line) of (a) the spin bare function at q = Q = (tr,7r) and 
(b) the 77 -mode bare function Xrj <1 =(0,0) at T = OA" in the Fermi 
liquid phase. A broadening S of 5meV was employed. 

references'*® and^^), we find; 

= r^Tr[G(/e-Fia),k + q)-(5(ie,k)] , (11) 

£,k 

where e{co) is a fermionic (bosonic) Matsubara frequency, T 
the temperature and Tr means Trace of the Green function ma¬ 
trix G (see expression (2)). G is the Green function matrix de¬ 
duced from the inversion of the matrix (2). From the figure 5 
a), and assuming that is frequency and momenta indepen¬ 
dent, the ladder diagram x^ can be determined (see"*® and^^) 
in the form: 

4(k,q,a),e) = 7’ ^ 
ffli-qi 

Tr [G(/e-l-i®-F/a)i,k + q + qi -fqc) • G{—ie — !®i,—k —qi)] 

( 12 ) 

where coi is a bosonic Matsubara frequency. Note that we put 
J-q outside of the expression (12). We do the internal Matsub¬ 
ara frequency and momentum summation by doing the change 
of variable k = k + qi and cb = coi + e. Then, Xq (k, q, to, e) 
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can be simplified as (fl, becomes : 

X°(q,®) = T ^Tr[G(ia) + ia),k + q + qc) -^(-/oij-k)] 


The bare polarization in the spin and Tj channels in the 
SC/CDW coexisting phase are presented in Fig.6 (with 
= 63 meV and = 35 meV) and in the Fermi liquid 

phase in Fig. 7. The bare polarization in the spin channel Xs 
develops a gap in the SC/CDW phase Fig.6 (a). The presence 
of a threshold in Imj® at the frequency CO « 60meV allows the 
emergence of a resonance of the collective response below the 
gap. This gap closes in the Fermi liquid phase as seen in Fig. 7 

(a) . 

On the other hand, the Tj-mode couples an electron with 
momentum k + qe to a counter propagating electron with mo¬ 
mentum —k. Consequently, Xq has an opposite sign of the 
one in the spin channel in both Fermi liquid and supercon¬ 
ducting state (see Fig.6 (b) and Fig.7 (d)). In the rj channel 
X^ develops a gap in the SC/CDW phase Fig.6 (b). This quasi¬ 
particle gap closes in the Fermi liquid phase as seen in Fig.7 

(b) . We clearly see a gap with a threshold occurring around 

(0 « 60meV which is the value of Note that a gapped 

quasiparticle continuum develops at low frequency in x^ as 
in Xs- The reason is that we evaluated x^ at q = 0 while the 
CDW order has a non-zero ordering vector qp. For compar¬ 
ison, Xs evaluated at finite momentum q = Q which also 
results in a gapped quasiparticle continuum. 


2 . 


Calculation of the term 


4 ’^ 



FIG. 8. (Color online) a) The Raman bare polarization bubble with 
the Raman vertex in the X symmetry, b) Xcoul Coulomb 

screening term resulting from the RPA series of diagram. Note that 
the dark triangle in the diagram denotes the presence of the Raman 
vertex. 


D. The full Raman response 


We now argue that the t]-collective mode is seen in the Aig 
Raman channel only. The study of the Raman susceptibility 
requires a careful examination of the symmetries of the sys¬ 
tem. These symmetries are taken into account by considering 
vertices in the photon-matter interaction different from unity 
(see Fig. 4 and 8). Three symmetries are typically considered 
for the Raman vertices in the cuprate superconductors which 
write within the effective mass approximation: 








; y®'.? = 


' d% 

_d^kx 


1 

2 dkxdky 


+ 


dkydkx 


d% - 

d^ky 


(15) 


The ladder diagram 4’^ can be determined (see^® and^^) 
in the form: 

4'‘(q,®)=-T^ E ■^h(®Fqi)44+q+q,X 

,k,qi 

+ kH-q) •G^^c(,£^k)x 


where A^g probes the whole Brillouin zone. Big the antinodal 
zone and Z? 2 g the nodal zone. Here we assume that there is no 
response in the A 2 ^ symmetry*® (■jAzg = q). In presence of the 
TJ collective mode, the full Raman susceptibility writes in the 
A symmetry: 


— -yt* -I- 4- 

ARaman AyXyX ' Acoul ' At] 


(16) 


qA:dw (jg -p , k + q + qi -f qc) X 

G'^^c (_/£ _ jgji _k - qj) 


(14) 


Considering our initial approximation, we can factorize by 
the interaction Jrj and we determine, from the Green func¬ 
tion matrix (2), that — /OTi, —k — qi) = 4- 

/a)i,k + qi) and G^™iv(;e _j_ j® k + q + qi Tq^) = 

G^cow,*(^i£ -|-/®-|-i®j,k + q + qi). We do the internal Mat- 
subara frequency and momentum summation by doing the 
change of variable k = k + qi and O) — (0i + £■ 

This susceptibility is proportional to the product of both SC 
and CDW order parameters Acdw-^sc- Consequently, the Tj- 
mode susceptibility 4 * only exists in a coexisting SC/CDW 
phase. A careful study of the symmetry of the Tj-mode re¬ 
sponse is done in the following section. 


where ’■^e bare Raman response (Fig. 8 a)), Xcoul 

Coulomb screening (Fig. 8 b)) and x^j the Tj-mode response 
(Fig. 4 a)). 

The value of 4 depends on the momentum dependence 
of the product y^G^coiv . q^sc , xhe product of the two 
Green functions • G'^-^c is proportional to the product 

Asc-Acdw which is proportional to the square of the d-wave 
factor. The sum over the internal momenta in the FBZ of the 
square of the tZ-wave factor is finite. The sum over the internal 
momenta in the FBZ of the Raman vertex 7 ^*® and 7 ^^® van¬ 
ishes {Y,k 7^= Lit = 0) while it is finite for the Raman 
vertex ® (£^ Yk ^ ^ The sum over the internal momenta 
of the product 7 ^G^°'^ • G^®^ therefore vanishes in the Big 
and B 2 g symmetry while it is non-zero in the Aig symmetry. 
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Consequently, the T 7 -mode response only exists in Aig sym¬ 
metry. 

In addition long-range Coulomb interaction U ^ l/q^ plays 
an important role in screening the Raman in the A channel . 
Doing the summation over all the coulomb diagram (Fig. 8 b), 
one can find easily: 




(17) 


where is the Raman susceptibility with one of the ver¬ 
tex put to unitiy and Xn the bare polarization bubble. In the 
limit g' —> 0, the "Coulomb screened" susceptibility simplifies 
as xLi = In the Ajg symmetry, the contri¬ 

bution of Xcoui cannot be neglected*^. Its contribution will 
screen partly the bare Raman susceptibility Con¬ 

sequently, the Raman response in the Aig symmetry writes 
Xplan = +xtaui + X^''- Around optimal doping, we 

will see that the contribution of the T 7 -mode (last term) is of 
the same order of magnitude as the bare screened A Raman 
response (first two terms). By increasing the doping, because 
of the weakening of the CDW order we expect a decrease of 
the Tj-mode intensity and the Aig spectra will be dominated 
by the bare screened Raman response. 

The Coulomb screening can be neglected in the Big and 
B 2 g channels for symmetry reasons*^. Consequently, the full 
Raman susceptibility is simply given by the unscreened re¬ 
sponse: xLman « with A = Big{B 2 g). 


III. RESULTS AND DISCUSSION 

The theoretical spectra are calculated without solving the 
self-consistent equation derived from the spin-fermion model. 
As a consequence, we need to adjust the amplitude of the or¬ 
der parameters as well as the value of J and V to reproduce 
the experimental data. The parameters values are adjusted in 
the following way : the value of A^^^ is chosen to reproduce 
the Raman peak frequency in the Big channel. The value of 
^CDW chosen to adjust the intensity of the resonance in the 
Aig channel. Note that the value of does not influence 
the frequency of the Tj-mode. We adjust the value of J to fit 
the neutron resonance frequency and the value of V to fit fre¬ 
quency of the Raman resonance in the Aig symmetry. A full 
self-consistent determination of the CDW and SC order pa¬ 
rameter amplitude is not done here as our goal is to identify 
the main theoretical features necessary to explain the Raman 
resonance in the Aig symmetry. 


A. Results 

The imaginary parts of the neutron spin susceptibility (spin 
response) and the full Raman susceptibility in the Aig sym¬ 
metry (Raman response) are presented in (a) and (b) panels of 



FIG. 9. (a) Calculated spin response, (b) Experimental (solid line) 
and calculated (dashed line) Raman response in the Aig symmetry 
and (c) in the Sig symmetry at optimal doping p = 0.16 (5 = WmeV, 
2Ag(- = 63 meV, = 45meV). In (b), the dotted line presents 

the calculated bare screened Raman response without the Tj mode 
contribution, Phonon lines have been sub¬ 

tracted out of the experimental Raman spectra for clarity in (b) but 
not in (c). The arrows in (c) indicate the location of the phonon lines 
superimposed on the electronic background. 


Fig. 9 for optimal doping (p=0.16) with 2Ag(g = 63meV and 
= 45mey. The theoretical curves are obtained for 
J = 107.5 meV and V = 54.5 meV. To the best of our knowl¬ 
edge, the amplitude of the CDW order parameter has 

not been evaluated yet in the experiments. However, in our 
effective model, the value of is reasonable compared 

with the SC order parameter amplitude. 

For this set of parameters, the calculated spin response ex¬ 
hibits a sharp peak at O) = 364cm^^ at the same energy than 
the energy of the calculated Aig Raman resonance. Note that 
the frequency dependence of the CDW order Acdw cuts the 
SC coherence peak at 2Ag(^ ~ (510a«^* (63mey) in the Aig 
channel but leaves the Big channel unaffected (see panel (c)). 
Fig. 9 (b) and (c) show a qualitative agreement between the 
peak energies in the Aig and Big calculated Raman spectra 
(red curves) and the experimental Raman spectra obtained 
from Bi2212 single crystals (black curves)^ \ We note that 
in our approximation, the calculated susceptibility vanishes in 














9 



wave number co (cm'^) 

FIG. 10. (a) Calculated spin response, (b) Experimental (solid 

line) and calculated (dashed line) Raman response in the symme¬ 
try and (c) in the B\g symmetry for the overdoped regime p = 0.23 
(5 = WmeV, 2Ajp = 25meV,2A^^^/ = OmeV). Phonon lines have 
been subtracted out of the experimental Aig Raman spectra in (b) for 
clarity. 


the metallic phase. This implies that we cannot reproduce the 
susceptibility at high energy for (O >2A®(^. More realistic elec¬ 
tron self-energies are needed in order to reproduce the spectra 
at high energy. 

In the over-doped regime (p=0.23), the value of 7 = 65 meV 
decreases and the value of V is the same than at optimal dop¬ 
ing. Experimentally above optimal doping, the frequency of 
the INS resonance at Q decreases with overdoping and scales 
with In absence of experimental data, we choose to ht the 
neutron scattering peak by assuming that the scaling between 
the neutron resonance frequency and Tc postulated in^ is still 
valid at hole doping p — 0.23. At this doping, the INS reso¬ 
nance frequency should be at 187 with Tc = 5QK. Note 
that persistent spin excitation have been observed in a huge 
range of hole doping by resonant inelastic X-Ray scattering 
experiments^^’^® which favours the existence of neutron reso¬ 
nance at Q vector at high doping. The maximal value of Asc 
and Acdw weakens and we take into account this decreasing 
by putting 2A®(^ = 25 meV and 2 A^q^ = OmeV. The decrease 
of A®^^ is due to the larger distance to the AF QCP at this 
doping as well as the proximity of a change of Fermi surface 


topology^ The vanishing of the CDW order in overdoped 
regime has been observed by STM^*. 

For these parameters the calculated spin and Ai^ Raman 
responses exhibit peaks at different energies (187cm^* and 
205 cm^* respectively) ((see Fig. 10 (a) and (b)). The perfect 
matching between Neutron and Raman peak energy is thus 
only verihed close to optimal doping level. The theoretical 
Raman response capture well the energies of both the Big SC 
coherence peak 2A®(^ and the Ai^ peak detected in the experi¬ 
mental spectra shown in Fig. 10 (b) and (c). 

The rj-mode also exists in the overdoped regime but the de¬ 
creasing of Asc and Acdw completely weakens its intensity. 
The 2A SC coherence peak, which is not detectable in the op¬ 
timal doped regime (see Fig.9 (b)), increases its intensity with 
overdoping and become detectable in the Aig channel around 
p = 0.22. Beyond this doping, the CDW order parameter am¬ 
plitude vanishes and the q-mode contribution to the Ai^, Ra¬ 
man spectrum vanishes. At this doping levels the Aig spec¬ 
tra are entirely given by the bare screened Raman response 
which peaks at a different energy than the neutron resonance. 
To summarize, the physical mechanism behind the Aig reso¬ 
nance peak is different below p = 0.22 than above. Below 
this doping, the Aig resonance originates the 77 -mode. Above 
this doping, the 77 -mode vanishes with the CDW order and the 
peak in the Aig channel is a 2A SC coherent peak located at 
the same energy as the one in the Big channel. This is qualita¬ 
tively consistent with Raman experiments, more neutron data 
in the overdoped regime would be desirable to conhrm this 
picture. Note that in the Fig.9 and 10, we adjust the units 
to superimpose the experimental and calculated data. Con¬ 
sequently, we cannot comment on the spectral weight of the 
results. 

Before proceeding with the discussion, we briefly com¬ 
ment on the choice of parameters. In the AF state, ab-initio 
calculations determine a superexchange coupling = 
125meV^^~^^. An experimental value of superexchange cou¬ 
pling J^p has also been deduced from Raman scattering 
experiments in Bi2212 compounds and evaluated at 7 « 
125mey®4.65 

in the AF phase. With hole doping, it is rea¬ 
sonable to assume that the superexchange coupling decreases 
7 < Jaf^^ because of the screening caused by the presence 
of holes. The values of 7 we deduce from the fits (see Fig.9 
and Fig. 10) are comparable with the calculated and the ex¬ 
perimental values, and are coherent with a decrease of the 
superexchange coupling with the hole doping. The value of 

V is weaker than the superexchange 7. To the best of our 
knowledge, this value has not been evaluated yet. The non 
zero value of the V parameter implies that the charge channel 
is important to explain the Aig Raman resonance. The fre¬ 
quency of the Ai^ Raman resonance depends on the values of 

V and 7 since it depends on 7^) —3J — V. Note that the Aig 
Raman frequency is less sensitive to the value of V than to the 
value of 7. Note that our effective model does not involve the 
strong correlation existing in the system. The value of 7 and 

V reported here are then effective parameters. 
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B. Discussion 

In this part, we discuss our theoretical results and the main 
experimental features that characterize the Raman spectra in 
the Aig and Big channels in the cuprates. 

1) The Raman experiments in the Big symmetry well estab¬ 

lished since a long time the presence of the superconducting 
coherent peak at the frequency of twice maximum of the d- 
wave SC gap Note that the value of can be 

determined by STM experiments^ ^ Our theoretical model re¬ 
produces well this superconducting coherent peak for a large 
range of doping (see Fig 9 and 10)^^. 

2) In Raman experiments at optimal doping, the Raman Aig 

resonance appears at a frequency below the superconducting 
coherent peak frequency 2 A^(-^“^^ and at an energy very close 
to the INS resonance frequency at Q = {n, Our model 

explain this resonance by a Tj-collective mode which rotates 
the SC ground state into a CDW state. The energy of this reso¬ 
nance depends on the interaction amplitude J and V as shown 
in Fig. 9 (a) and (b). Althought our model can reproduce the 
nearly identical energies of the Aig Raman resonance and the 
neutron resonance, the equality is obtained by adjusting the 
parameter V. Nevertheless our model is able to reproduce this 
equality for a reasonable set of parameters 

3) In the overdoped regime, the Raman Ajg resonance is ex¬ 
perimentally observed at an energy very close to the supercon¬ 
ducting coherent peak frequency 2A®(^. In our calculations, 
the 77 collective mode vanishes in this regime because of the 
disappearance of the CDW order. Consequently, only the bare 
screened susceptibility contributes to the Raman response and 
the resonance becomes a superconducting coherence peak 
as shown in Fig. 10 (a) and (b). 

4) In the underdoped regime, the transition temperature of 
the CDW order is greater than the superconducting critical 
temperature. Hence we can consider that the CDW order pa¬ 
rameter amplitude becomes larger than the SC one A^^j^ > 
Aj(^. In our theory this regime implies a decreasing of the 
T 7 -mode resonance intensity which will become undetectable 
at low doping. Existing Raman data seems to be consistent 
with a disappearance of the Aig in the underdoped regime^^. 
The description of the underdoped regime however requires a 
description of the pseudogap phase that we did not take into 
account in the present work. Deeper considerations about the 
Raman scattering spectra in the underdoped regime will be 
considered in details in a forthcoming work. 

5) To fit the Raman resonance in the Aig symmetry, we pro¬ 
pose an explanation based on the coexistence between two 
orders with a d-wave symmetry. This d-wave symmetry re¬ 
mains crucial to explain the resonance in the Aig symmetry 
only. Moreover, the d-wave form factor of both SC and CDW 
orders has been observed by STM^^. To explain the small am¬ 
plitude of the 2 A superconducting coherent peak in the Aig 
symmetry, it is necessary to include a frequency dependence 
of the order parameters that can be predicted by the spin- 
fermion model. This frequency dependence suppresses the 
superconducting coherent peak from the Tj-mode response. In 
the bare Raman response, this superconducting coherent peak 


is strongly screened by the long-range Coulombian effects. 
This screening allows the observation of the Tj-mode in this 
channel. 

6 ) Experimentally, the Aig Raman resonance has been ob¬ 
served in the SC state and disappear at . The Tj-mode 
calculated here depends on the magnitude of the CDW and 
SC order parameters and then we expect to see its disappear¬ 
ance for temperature below Tc in the overdoped part of the 
phase diagram (where A®^ > This is a limitation of 

the approach developed here. However, we argue that the Tj- 
mode exists between SC order and all type of charge orders. In 
our simplified model, we neglect the contributions that could 
come from the PG phase. In the framework of the interpre¬ 
tation developed in Ref.®®, where the PG phase is interpreted 
as a charge order, one could observe the emergence of such rj 
resonance at higher temperature than the CDW critical tem¬ 
perature. 


IV. CONCLUSION 

To conclude, we propose a coherent scenario to explain 
the resonance peak in the Aig symmetry seen in Raman scat¬ 
tering. This scenario is based on the existence of coexist¬ 
ing ^f-wave SC and (f-wave CDW orders in the SC phase 
of the Cuprates. The ubiquitous proximity of the CDW and 
SC states in the underdoped phase of the cuprates generates 
a collective “rj-mode” allowing rotations between those two 
states. The J] collective mode results from the coupling be¬ 
tween the d-wave SC and the d-wave CDW state. This mode 
produces a response in Raman scattering spectroscopy solely 
in the Aig symmetry and matches the spin triplet resonance 
at q==Q=(;r,;r) at optimal doping. Since the Tj-mode is 
a charge ±2 spin zero spin singlet excitation other probes 
like Electron Energy Loss Spectroscopy and Resonant X-Ray 
techniques are also likely to show the resonance. 
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Appendix A: Supplementary Informations about the 
spin-fermion model 


out the bosonic degrees of freedom, we neglect the spinwaves 
interaction (m = 0). The partition function then writes 


In this section, we present the calculation of the effective 
Lagrangian derived from the initial spin-fermion Lagrangian. 
The goal of this section is to give more physical details for the 
completeness of the paper. 


1. The spin-fermion model and the effective boson propagator 

The starting point of the theoretical model is the spin- 
fermion model^®’^®. Note that the models presented 
simplify the Fermi surface (FS) in Eight Hot Spots (EHS) (the 
hot spots are points of the Eermi surface separated by the vec¬ 
tor Q = {n, n)). In the following, we develop the model intro¬ 
duced in^® in the whole Brillouin zone (BZ). 

In the model^®, the physics of cuprates results from the cou¬ 
pling of the fermions with a bosonic spin mode which gener¬ 
ates at the Antiferromagnetic(AE)/PseudoGap(PG) Quantum 
Critical Point (QCP) a pseudogap in the fermionic dispersion. 
This new state, different from the AE one, is a superposition 
of a ^f-wave superconducting (SC) state and a Charge Den¬ 
sity Wave (CDW) order. This bosonic spin mode (also called 
paramagnons) survives out of the AE phase and become the 
new glue between electrons that leads to the PG and super¬ 
conducting (SC) phase. To describe the interaction between 
this bosonic spin mode and fermions, we start from the La¬ 
grangian L = Lc + L^, where 

Lf = c* ((9 t- f -|-Z0(7)c (Al) 

L(j) = — + m)(j) + (A2) 

2 vj 2 

The Lagrangian L describes the coupling of the electron (rep¬ 
resented by the fermionic field c) to spin waves (represented 
by the bosonic field 0). The dynamic of the fermions and their 
coupling with the bosonic mode is described in the Lagrangian 
Lc- The strength of this coupling is tune by the parameter A 
which couple c and 0. The fermionic dispersion is described 
by 1 *^ = e^. — /r where is the fermionic dispersion and jj. the 
chemical potential and afford us to describes the Eermi surface 
of the system. 

The bosonic spin mode physics is described in the La¬ 
grangian L^. Vi is the velocity of the spin waves, (O the 
bosonic Matsubara frequency, the paramagnon mass m that 
describes the distance to the QCP (m > 0 in the metallic side 
and vanishes at the QCP) and Q = {7t,n) the ordering wave 
vector in the AE phase. 

The spin wave propagator -f (q —Q)^ 

is renormalized by the particle-hole bubble (see details in 
Ref ^^). This renormalization leads to consider the effective 
spin wave propagator D^ff that writes j^ = (tI®I + |kp + 
m) where y a phenomenological coupling constant and k has 
formally been shifted by Q. 

In the following, we simplify the notation as k = {ico,k), 
where ico are fermionic Matsubara frequencies. To integrate 


Z = J d[^>]expi-So-Sint) (A3) 

= (A4) 

k,G 

Sint = — 52 ■^‘l'^ta'^k+Q+q,acl, (A5) 

k,k' ,q,a 


where the bare propagator is 


^Ok 


( ico - 
0 
0 

V 0 


0 


l(0~ 


o-k+qc 


0 

0 



(A6) 


and the spinor field '¥k = ^). More¬ 
over, 7^' = 4D^'/3A^, a = labels the spin and Qc 

stands for the charge density wave ordering vector. 


2. The mean field decoupling and the effective Lagrangian 


After the integration on the bosonic degrees of freedom, 
we decouple the two-body action Sint using mean field decou¬ 
pling. Thus, we decouple the interaction Lagrangian in the 
SC channel and in the charge channel by introducing the two 
order parameters A*, and AcDW,k- Ak is the superconducting or¬ 
der parameter which describes the 7-wave SC order. AcDW,k is 
the CDW order parameter which describes the 7-wave charge 
order with ordering vector Qc. Note that in the reference^®, 
these two order parameters are related to the SU(2) symme¬ 
try which is verified at the hot-spot. The SU(2) symmetry 
between SC and CDW order exists if their energy level is the 
same (the energy splitting Ak— Acnw.k vanishes). In the whole 
Brillouin zone, the SU(2) symmetry is verified around the hot¬ 
spots (at low paramagnon mass) and at the zone edge (at high 
paramagnon mass)^^. 

Applying the Hubbard-Stratanovitch transformation to the 
interaction action Sint, the partition function becomes 


/ d['¥]d[A,AcDw]exp{-So Smt,e ff) 
J d[A,AcDw]exp{-Stj) 

~ ~ ^ l-^q ^CDW,k^CDW,k+q + Jq 
k,q,o 

Sint,eff = St,-Y,^laMk^k,a, 
k,o 


{Al) 

(A8) 

(A9) 


with k = k + Q and the matrix A4 is 
ihk 


Mk = 




,mk = 


—AcDW,k —Ak 

~\+qc ^CDW-k 


(AlO) 


Einally, the effective Lagrangian we will use in the follow¬ 
ing will be defined as = Gq^* — Mk- One can integrate 
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used from respectively a diode pump solid state laser. Mea¬ 
surements in the SC state have been performed at lO/T us¬ 
ing an ARS closed-cycle He cryostat. Raman study was per¬ 
formed on Bi 2212 single crystals with two distinct levels of 
doping p = 0.16 {Tc = 9QK) and p — 0.23 {Tc = 52K). The 
level of doping was controlled only by oxygen insertion (see 
Ref 51). TheBi^ channel was obtained from crossed polariza¬ 
tions at 45*^ from the Cu-O bond direction. The Aig channel 
was obtained from parallel polarizations along the Cu-O bond 
(given Aig-t Big) and normalized subtraction of the Big chan¬ 
nel. All the spectra have been corrected for the Bose factor 
and the instrumental spectral response. They are thus propor¬ 
tional to the imaginary part of the Raman response function 
Imx^(®). (X =Aig or Big) 


FIG. 11. The frequency dependence of the superconducting order 
parameter Aj^ at the hot-spot calculated from the mean field equation 
(All). The frequency form factor which is a Gaussian with a half¬ 
width of Aq which is the value of the SC order parameter at zero 
frequency. 


over the fermionic degrees of freedom and find the mean field 
equations as in Ref.^^: 

M(o.k = T ^/k-k'-gGm'k'- (All) 

ffl'.k' 

The matrix equation can be projected onto the different com¬ 
ponents. The exact solution of these equations goes beyond 
the goal of the paper. However, we directly see that the order 
parameters depend on the frequency as well as on the momen¬ 
tum. 


3. Frequency dependence at the hot-spot 

The frequency dependence is a key player of the theory. 
Then we solve the equation (All) for the superconducting 
component A® t only, putting CDW order parameter to zero 
(A^^ = 0). Considering a temperature of T = 0.001 and a 
paramagnon mass to m = 10^® with 7 = 10^^ and X = 12, 
the frequency dependence of the SC order parameter at the 
hotspot is shown in Fig. 11 . 

As shown in Fig. 11, the Gaussian frequency form fac¬ 
tor can reproduce the frequency dependence of the SC order 
parameter. In the main text, we assume that the frequency de¬ 
pendence of the CDW order parameter can be fit by the same 
form factor. The calculation of the frequency dependence of 
the order parameters at each doping and each temperature is 
left for later work. 


Appendix B: Raman experiments 

Raman experiments were carried out using a triple grating 
spectrometer (JY-T64000) equipped with a liquid-nitrogen- 
cooled CCD detector. The 532 nm laser excitation line was 
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